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1. INTRODUCTION 

Brushless dc motor has windings in the stator and permanent magnet in the rotor. As it does not 
have brushes and commutator the efficiency of the motor increases with decrease in the ohmic loses. 
Brushless dc motor are available in single phase, two and three phase configuration. But three phase 
configuration are commonly used. Generally, Hall Effect sensors are used to sense the position of the rotor. 
But in this project the rotor position is determined from excitation current and voltage across the stator 
windings. The BLDC motor requires inverter as there is no brush and commutation arrangement. 
The inverter uses transistor for low drive application and Thyristor for high power drives. MOSFET/DIODE 
switches are used in the inverter. The commutation sequence to the MOSFET/DIODE depends on the rotor 
position. BLDC motor require less maintenance, generate less noise, have low inertia, have high efficiency, 
has long operating life. 

The controller is used to improve the dynamic response of the system. The PID and FOPID 
controller parameters are tuned by genetic algorithm [1], more optimization techniques are used to tune 
the coefficients ofFOPID and PID controller [2]. In this paper the process to control is speed of the BLDC 
motor. The BLDC motor is modeled [3, 4]. The FOPID controller transfer function is formed by using 
FOMCON toolbox [5, 6]. The genetic algorithm optimization to tune PID and FOPID is done by suing 
optimization toolbox. 
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2. SPEED CONTROL OF BLDC MOTOR 

The block diagram for speed control of BLDC motor is shown in Figure 1. The rotor position and 
speed of the BLDC motor is computed from the voltage and current supplied to the stator windings of BLDC 
motor. Each phase winding is excited in a sequence to run the BLDC motor. The phase winding of the BLDC 
motor is excited in a sequence by sequence commutation of the MOSFET/DIODE in the inverter [7, 8]. 








reference 
speed 


















VARIABLE 
DC 
SOURCE 





THREE 
PHASE 
INVERTER 


FOPID OR 
PID 
CONTROLER 














position 
ofthe 
rotor 


current 
speed 





POSITION 8 TO 
GATE SIGNAL 


Figure 1. Block diagram of the proposed system 


That switching sequence is given to the inverter from position to gate signal block as shown in 
Figure 1. The speed of the BLDC motor is directly proportional to the voltage. Therefore, to control 
the speed of the BLDC motor, the input voltage must be controlled. The speed computed is compared with 
the reference speed by the comparator and the error signal is obtained. The error signal is given as input to 
the controller. The output of the controller is control signal. This control signal is given as input to 
the variable DC source. There by the input DC voltage is controlled based on the BLDC motor speed. 


2.1. PID controller: 

The PID controller consists of the proportional, integral and derivative controller. He proportional 
controller gives proportional response of the input error value. The integral controller is proportional to 
magnitude and the duration of the error. It improves the settling point and eliminates steady state error. 
The integral term gives response which also depends on the sum of previous error values. The sum of 
the previous error value is multiplied with the gain ki. The derivative control gives the slope of the error for 
specified duration and multiples it with the gain kd. The sum of this proportional, integral and derivative 
control gives the PID controller. The control signal of PID controller in time domain. 


u(t) =k, e(t) + ki feoai + kg Ten) 
0 


Taking Laplace transform of above equation, the transfer function of PID controller is given by, 


U(S) _ 


ki 
E(S) ky + P kas 


Ge (S) = 
Where, kp is proportional gain, ki is integral gain, kg is derivative gain, e(t) is an error signal. 
The block diagram of PID controller is shown in Figure 2. The parameter kp, kj,kg have to be tuned 


by using genetic algorithm. The PID controller can implement in MATLAB/Simulink by PID block available 
in continuous time toolbox. 
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Figure 2. Block diagram of PID controller 
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2.2. FOPID controller 

Fractional order PID controller is commonly in the form PI*D". In this the integrator and 
differentiator is in the order of à and u respectively [9]. The time domain representation of FOPID control 
signal is, 


u(t) = kpe(t) + k,D~*e(t) + kgD"e(t) 


Where, D~* - Fractional order integrator, D} - Fractional order differentiator. Taking Laplace transform for 
above equation, the transfer function of FOPID controller is given by, 


U(S) _ ki 


G.(S) = E(S) | ky + så 


+ kas" 


Where, A - Order of integrator, u - Order of differentiator. The block diagram of FOPID controller is shown 
in Figure 3. 
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Figure 3. Block diagram of FOPID controller 


3. DESIGN OF BLDC MOTOR 

The rotor and shaft of the BLDC motor is assumed to be rigid [10]. The mode is assumed as having 
viscous friction model. Therefore the friction torque is proportional to the shaft angular velocity. Based on 
Newton second law and Kirchhoff s voltage law, three phase star connected BLDC motor can be described 
by following (1-4). 


d Om 





Te = bOm +J +T, (1) 
Vap = Rig — ip) +L” + e, — (2) 
Voc = R(ip — ic) +L? + e,- e, (3) 
Voa = Ric — ig) +L + e,- eg (4) 


Where Te,Om,; b, J, TL, V, R, L, I, e are electrical torque, mechanical rotational speed, viscous friction 
constant, the rotor inertia, mechanical load torque, phase to phase voltage, phase resistance, phase current, 
phase inductance and phase back emf respectively. 

The voltage and current equation is given by (5-6), 


Vap + Vic + Vea = O (5) 
lap + Ipe + lca = O (6) 


To simplify the modeling, only two voltage equation are need, they are (7-8). 


dia 


2V ab + Vpc = 3Rig + 3L dt 


26, Se), =e. (7) 
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-Vap + Voc = 3Rip + 3L% + 2ep — ea — €, (8) 
The torque generated by the BLDC motor is given by (9), 
; A ; dOm 
Te a (eala + eplp + ecte M (9) 


The trapezoidal back emf can be written as (10-12) : 


Ke dðm 

€a = > [ra(ðe) (10) 
Ke dðm 2 

Cp = z ay Tae = 5 (11) 
Ke dðm 4 

ec = 7 [rale E =) (12) 


Wherek,, Oe is the back emf constant and electrical angle respectively. 
Oe = pO» 
Tra(@,) is the trapezoidal waveform unction and one period of the function can be described as follow (13): 
2m 
1, 0 < 0; < = 
6 270 270 
Tra.) =4 1-2(0,-=), <0, <n (13) 
6 57T 57T 
-1+°(0,-=), <0, < 2r 
Substituting (10-12) in (9) we torque as follows (14): 
Ke , 2n\ . 4T) . 
Te = “£[Tra(®,)ig + Tra (0. — 7) ip + Tra(®, - =) ic] (14) 


The angle of the rotor and speed of the rotor is related as follows (15-16): 


dOm 
Wap (15) 


We = PWm (16) 


Where, w, - Electrical speed, w,, - Mechanical speed. From the equation (1) and (15), we get the speed 
equation as (17): 


Wm = — (17) 

From (7), The current i, is given as (18): 
_ 2VaptVpc—2eqateptec—3Rig 

Mm 3Ls (1e) 
From (8), as 

. _ VaptVpc—2eateptec—3Rig 

— —— — (19) 
From (6), The current i, is given as (20): 

i = ig — İp (=o) 
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4. INVERTER TOPOLOGY 

The inverter used here is a three phase inverter. It consist of three leg and six MOSFET/DIODE. 
It is implemented by using universal bridge available in toolbox. The gate signal in this universal bridge 
decides the MOSFET commutation. The inverter is shown in Figure 4 And for Circuit diagram of dc motor is 
shown in Figure 5. 
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Figure 4. Three phase MOSFET/DIODE inverter Figure 5. Circuit diagram of DC motor 


5. SYSTEM MODELLING 
Transfer function of BLDC motor: Using Kirchhoff’s voltage law (21), 


y=Ri+L<+e (21) 


Where V,, R, I, L, e are the source voltage, resistance in ohm, current to the motor, inductance of the coil, 
back emf respectively. The torque of the dc motor from second law of motion is given by (22), 


Te = kW +J Z= +T, (22) 


Where Te,kf,J,Wm.T, are the electrical torque, the friction constant, the rotor inertia the angular velocity and 
the supposed mechanical load respectively. The electrical torque and the back emf could be written as, 


e=k,w,, and Te = ktWm 


Where k, is back emf constant and k; is torque constant. From equation (21) and (22). The transfer function 
for dc motor is given as (23): 


_ Wm _ kt 
G(s) = Vs  S?JL+(RJ+kfL)s+kfR+kekt (23) 
Considering the following assumption, 
— The friction constant is small (k; = 0) 
- RJ >K-L 
kekt > ke R 


And multiplying by — xiin (23). The transfer function of dc motor is given as (24): 
Kekt R 


1 


— ke 
G(s) = TmTeS2+tmst1 (24) 
: i RJ 
Where mechanical time constant Tm = ET 
eñt 


. š L 
Electrical time constant T, = = 


But for the BLDC motor 
_y RL 
Tm = Ly Kt 
L 
Te = R 
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Therefore, since there is a symmetrical arrangement and a three phase the mechanical and electrical constant 
become (25), 





_ 3RJ 
m keki 
L 
Te = 7r (25) 


Considering the phase effect (26), 


3RoJ 
_ p 
Tm = ker-i 
V3 





(26) 
kt 


k 
Where, k, = a From the MAXON EC 45 FLAT BLDC motor datasheet, the value of phase resistance, 


phase inductance, rotor inertia and mechanical time constant are shown in Table 1. 


Table 1. MAXON EC 45 FLAT BLDC motor datasheet 


Parameter Symbol Value Unit 
Phase resistance Ry 1.2 Ohm 
Phase inductance L 0.560 Mh 

Rotor inertia J 92.5 x 107° gcm? 
Mechanical time Tr 17.1 Ms 
constant 


Substituting the above values in (25), 


Ta = 155.56 X10 


Vv—-Secs 


From (26), ke = 0.076 





9 


rad 


Substituting the values of Tm, Te, Ke in (24) we get the transfer function of BLDC motor as (27): 





NC) = pace ea ep 
Transfer function of inverter: The transfer function of the inverter is given as (28): 

G(s) == (28) 
Delay time T;,, = 0.004 sec 

Gainkin = aes = 0.65 x" = 6,36 
Therefore (29), 
FO sa T 29) 
Finally, the transfer function of the system is given by (30), 

G(s) = G,(s) x G(s) 

Ee 83.38 (30) 


1.064107 9853 +7.106X10795s*+0.0211s+1 
The parameters are tuned for this system. 
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6. GENETIC ALGORITHM 


Genetic algorithm is an optimization technique. This technique is used to tune the parameter of PID 
and FOPID controller [11, 12]. The flowchart of the genetic algorithm is shown in Figure 6. 
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Figure 6. Flowchart of genetic algorithm 


It is an iterative process with the population in each iteration called a generation. First the initial 
population is generated. The population may range from hundred to thousand. From the population 
the fitness solution to each individual is chosen by objective function. In the selection process the portion of 
the existing population is selected. In the crossover/mutation, a new generation is generated. Finally, 
it checks whether it meet the criteria. If it meets the criteria it gets terminate otherwise the iteration repeats. 


6.1. Optimization of the PID and FOPID using genetic algorithm 
The genetic algorithm can be implemented by using optimization toolbox. In FOPID tuning 
the parameter k,,,k;,kg, 4, and u should be tuned to the proposed system [13]. The fitness function should be 


called as @(x) fun name(x). The lower and upper limit of the cost function is shown in the Table 2. 
The parameter of genetic algorithm is chosen as in Table 3. 


Table 2. Limit of FOPID parameter 


FOPID parameter Lower limit Upper limit 
Kp 0 2 
K; 0 2 
A 0 1 
Ka 0 2 
u 0 1 


Table 3. Criteria of genetic algorithm 


Population 200 
Initial population 10 
Selection Stochastic uniform 
Mutation Constraint dependent 
Crossover Constraint dependent 


In PID tuning, kp,ki,ka are the cost function. The limit of PID controller is chosen as in Table 4. 


Similarly by choosing criteria in optimization toolbox the parameter are tunes as shown in Table 5. 
The graphical representation is given as shown in Figure 7. 
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Table 4. Limit of PID controller 


PID parameters Lower limit Upper limit 
Kp 0 2 
K; 0 2 
Ka 0 2 


Table 5. Tuned value of parameter of controller 
Kp K, K; A u 
PID 1.98 1.782 1.869 1 1 
FOPID 1.904 1.7818 1.7061 0.1835 0.2217 








Figure 7. Comparison plot 


The graph generated while tuning the parameter of FOPID controller is shown in Figure 8 and for PID 
controller is shown in Figure 9. 
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Figure 8. FOPID step response by genetic algorithm Figure 9. PID step response by genetic algorithm 


7. SIMULATION AND RESULT 

The simulation model of the speed control of BLDC motor is shown in Figure 10. The PID 
controller is implemented by using PID block available in MATLAB. The transfer function of FOPID is 
generated by using FOMCON toolbox. The output speed response for PID controller is shown in Figure 10. 

An external disturbance is applied to the BLDC motor at 0.5 second. The PID controller controls 
the input voltage and the BLDC motor reaches the steady state at 8.5 second. Thus the settling time of 
the BLDC motor using PID controller is 8 second. The output speed response of the system due to FOPID 
controller is shown in Figure 11. 

An external disturbance is applied to the BLDC motor at 0.5 second. The FOPID controller controls 
the input voltage and the BLDC motor reaches the steady state at 1.3 second. Thus the settling time of 
the BLDC motor using FOPID controller is 0.8 second. Therefore the settling time of the speed due to 
the FOPID controller is 10.5 times faster than the PID controller. Comparison of PID and FOPID speed 
response is shown in Figure 12. 


Optimization and comparative analysis of PID and FOPID controller for BLDC motor (P. Vimala) 


182 o ISSN:2089-4856 


$ 


Figure 10. Speed response due to PID controller Figure 11. Speed response due to FOPID controller 





Figure 12. Comparison of PID and FOPID speed response 


In this section, it is explained the results of research and at the same time is given 
the comprehensive discussion. Results can be presented in Figures, graphs, Tables and others that make 
the reader understand easily [14, 4]. The discussion can be made in several sub-chapters. 


8. CONCLUSION 

The FOPID controller enhances the dynamic response of the BLDC motor compared to 
the conventional PID controller. The settling time of the speed response is very fast when we use FOPID 
controller. So the performance of the FOPID based speed controller is better and efficient. Stability and 
robustness of the BLDC motor can be improved. So this proposed controller can be extended to other 
industrial control system. 
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